Toward understanding the microscopic origin of nuclear clustering 



<N 

O 
<N 



J. Okolowicz,^ W. Nazarewicz,^'^'^ and M. Ptoszajczak^ 

^Institute of Nuclear Physics, Radzikowskiego 152, PL-31342 Krakow, Poland 
^Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996, USA 

^Physics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
^Institute of Theoretical Physics, University of Warsaw, ul. Hoza 69, PL-00-681 Warsaw, Poland 
^ Grand Accelerateur National dTons Lourds (GANIL), 
CEA/DSM - CNRS/IN2P3, BP 55027, F-14076 Gaen Gedex, France 

Open Quantum System (OQS) description of a many-body system involves interaction of Shell 
Model (SM) states through the particle continuum. In realistic nuclear applications, this interaction 
may lead to collective phenomena in the ensemble of SM states. We claim that the nuclear clustering 
is an emergent, near-threshold phenomenon, which cannot be elucidated within the Closed Quantum 
System (CQS) framework. We approach this problem by investigating the near-threshold behavior 
of Exceptional Points (EPs) in the realistic Continuum Shell Model (CSM). The consequences for 
the alpha-clustering phenomenon are discussed. 
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I. INTRODUCTION 

The immense richness of the nuclear many-body prob- 
lem stems from its genuine multi-scale character and un- 
derlying effective interactions that are strongly mediated 
by the nuclear medium. Further complexity is added by 
the OQS nature of the atomic nucleus, which requires 
a treatment of bound states, resonances, and the con- 
tinuum of scattering states within a unified framework. 
In standard applications of the nuclear SM, which is a 
standard tool for nuclear structure studies, the nucleus is 
described as the CQS where nucleons occupying bound 
orbits are isolated from the environment of scattering 
states 

To formulate the SM for OQSs, two frameworks have 
been proposed. The first one, the Gamow Shell Model 
(GSM)(2|, is the complex-energy CSM based on the 
Berggren ensemble [4]. The GSM, which is conveniently 
formulated in the Rigged Hilbert Space [5], offers a fully 
symmetric treatment of bound, resonance, and scat- 
tering states. The second one, the real-energy CSM 
in the Hilbert space [6- 9], is based on the projection 
formalism[lQ. .11]. A recent realization of this approach, 
the Shell Model Embedded in the Continuum (SMEC)[3, 
provides a unified description of structure and reactions 
with up to two nucleons in the scattering continuum us- 
ing realistic SM interactions. The CSM describes many- 
body states of an OQS using the effective non-Hermitian 
Hamiltonian in which the coupling to decay channels 
results in the anti-Hermitian component. The config- 
uration mixing due to the competition between Her- 
mitian and anti-Hermitian terms is a source of collec- 
tive features such as, e.g., the resonance trapping[12- 
[l4] and super-radiance phenomenon [15[ [l6j, m ultichan- 
nel coupling effects in reaction cross-sections [I3 and shell 
occupancies [iS^, modification of spectral fluctuations flQ^, 
and deviations from Porter-Thomas resonance widths 
distribution [l3, [2Q| . We shall argue that the appearance 
of cluster states close to the corresponding cluster emis- 



sion threshold is yet another consequence of continuum 
coupling [21]. 

Nuclear clustering is arguably one of the most mysteri- 
ous nuclear phenomena. The comprehensive understand- 
ing of its universal occurrence and properties is absent in 
the CQS formulation of the nuclear man y- b ody problem. 
The commonly used Cluster Model j22|, l23| is an a pos- 
teriori approach that assumes effective building blocks 
(clusters). A priori approaches, like the nuclear SM ap- 
proach, simply fail to predict cluster states at observed 
low excitation energies around cluster-decay thresholds. 
In this context, it is worth mentioning that there exists 
a formal correspondence between Slater determinants of 
single-particle harmonic oscillator wave functions - used 
in many CQS SM frameworks - and cluster- model wave 
functions [24]. However, this relation, explicitly utilizing 
the SU(3) dynamical symmetry of the harmonic oscilla- 
tor, is of limited use when it comes to quantitative de- 
scription of cluster states that requires a proper treat- 
ment of decay thresholds and asymptotic behavior of 
wave functions. In the most advanced approaches to 
cluster decay, the SM wave functions must be supple- 
mented with a cluster component to provide a quantita- 
tive agreement with experiment [25]. This failure of CQS 
approaches to describe cluster states is the central prob- 
lem in nuclear theory. 

In this study, we shall discuss collective features of SM 
eigenstates in the neighborhood of the one-proton decay 
threshold by analyzing realistic SMEC wave functions. It 
is expected that genuine features of the emergent collec- 
tive behavior of SM states are universal, regardless of the 
nature of the charged particle (proton, deuteron, alpha- 
particle, etc.) emission threshold, in analogy with the 
universal behavior of reaction cross sectionsjlTl, [26| and 
spectroscopic quantities [18]. In this sense, the collective 
behavior of SM eigenstates close to the one-proton decay 
threshold provides the insight into the formation mecha- 
nism of more complicated charged cluster configurations, 
like alpha-cluster configurations. The SMEC formulation 
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is particularly convenient for our purpose as it allows us 
to study the collective features of OQSs in terms of the 
competition of Hermitian (internal) and anti-Hermitian 
(external, through the continuum) mixing of SM eigen- 
states. 

In this work, collective features of CSM (SMEC) 
states are discussed by analyzing properties of excep- 
tional threads (ETs) of the complex-extended CSM 
Hamiltonian[27], i.e., the trajectories of EPs in the space 
of system energy E and continuum-coupling interaction 
parameter Vb[28]. EPs, corresponding to the entangle- 
ment of two sheets of eigenvalues by the square root sin- 
gularity, are results of an interplay between opposite ef- 
fects of Hermitian and anti-Hermitian parts of the CSM 
Hamiltonian[8]. The (£^, l^)-dependence of EPs around 
the proton decay threshold is studied to understand the 
near-threshold collective rearrangement of SM eigenfunc- 
tions, which is at the roots of clustering phenomenon. 

In Sec. [IT] we shall examine physical arguments that 
necessitate the use of OQS formalism in the description 
of nuclear clustering and recall its universal features. A 
brief summary of the real-energy CSM and the EPs is pre- 
sented in Sec. IIIII The genuine near-threshold properties 
of ETs are discussed in Sees. IIVI and |Vl using illustra- 
tive examples of SMEC eigenfunctions of ^^Ne and ^^S 
that are mixed through the continuum coupling to the 
£ = proton decay channel. In particular, we study the 
role of the so-called aligned state in which the collective 
strength is concentrated. In Sec. we investigate the 
dependence of the continuum-coupling correlation energy 
and the position of its maximum on fragment charges and 
the nuclear radius. These results provide an insight into 
collective effects in more complex near-threshold situa- 
tions involving decays of heavier charged particles such 
as ^H, ^He, ^He, and ^Be. Finally, the main results of 
this work are summarized in Sec. IVI[ 



II. CONFIGURATION MIXING NEAR THE 
DISSOCIATION THRESHOLD 

The cluster states are closely related to the nature 
of a nearby cluster-decay threshold. Ikeda et al. [29| 
remarked that alpha-cluster states can be found in the 
proximity of alpha-particle decay thresholds. Curiously, 
this insight has not raised much interest in the role of 
continuum coupling in generating correlations associated 
with effective cluster degrees of freedom. 

Actually, the conjecture of Ikeda et al. can be for- 
mulated more generally [21]; namely, the coupling to 
a nearby particle/cluster decay channel induces parti- 
cle/cluster correlations in CSM wave functions which are 
the imprint of this channel. In other words, the clus- 
tering is the generic near-threshold phenomenon in OQS 
which does not originate from any particular property of 
the nuclear forces or any dynamical symmetry of the nu- 
clear many-body problem. This generalized conjecture 
holds for all kinds of clustering, including unstable clus- 
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FIG. 1: The scheme of couplings in the multichannel rep- 
resentation of ^^C. The structure of a near- threshold Hoyle 
resonance in ^^C exhibits a strong imprint of the nearby 
^^C^^Be+a decay threshold. With the increasing excitation 
energy, subsequent decay channels open up, leading to an in- 
tricate multichannel network of couplings. 



ters such as dineutron or ^Be. Figure [T] shows an example 
of ^^C for which the lowest-energy threshold corresponds 
to the emission of an alpha-unbound ^Be. 

What can be said about the properties of many-body 
states in the narrow range of energies around the re- 
action threshold? Are they universal, independent of 
any particular realization of the Hamiltonian? The clus- 
tering phenomenon can, in fact, be traced back to ba- 
sic properties of the scattering matrix in a multichannel 
network [30j. The decay threshold is a branching point 
of the scattering matrix. For energies below the lowest 
particle-decay threshold at energy Ei , one finds the ana- 
lytic phase with a single solution which is regular in the 
entire space. At energies above Ei and below the next 
decay channel threshold at energy £^2, the new analytic 
phase corresponds to two regular solutions of the scat- 
tering problem, and so on. As a result, one obtains a set 
of analytic phases, each one with a different number of 
regular scattering solutions. These phases are separated 
by decay thresholds and, together, form the multichannel 
OQS. The coupling of different SM eigenfunctions to the 
same decay channel modifies regular scattering solutions 
in the proximity of each branching point. 

Such a network is illustrated schematically in the r.h.s. 
of Fig. [TJ For energies E below the first decay threshold 
Ei2c^^^^ the regular phase consists of elastic scatter- 
ing solutions. For E > Ei2q^^^ and below the second 
threshold at energy £^12^^11^^^, the regular phase con- 
tains two solutions coupled by unitarity, and the story 
repeats at higher excitation energies as more and more 
channels open up. The structure of any eigenfunction 
in this multichannel OQS depends on the total energy 
of the system which in turn determines the environment 
of decay channels accessible for couplings. With increas- 
ing excitation energy, the multichannel OQS network of 
couplings becomes more and more complex. The flux 
conservation (unitarity) implies that the mixing of SM 
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eigenfunctions changes if a new channel becomes avail- 
able. In this sense, the OQS eigenfunctions are not im- 
mutable characteristics of the system but vary with in- 
creasing energy of the system [8|. 



III. CONTINUUM SHELL MODEL 

The real-energy CSM formulation using Feshbach pro- 
jection formalism[10, 11] offers a convenient framework 
for a discussion of the competition between Hermitian 
and ant i- Hermitian components of the effective Hamilto- 
nian. The Hermitian term causes the repulsion of levels 
in energy and attraction of widths, whereas the anti- 
Hermitian term leads to energy-attraction and width- 
repulsion effects [8]. The resulting interplay between Her- 
mitian and anti-Hermitian couplings yields a complicated 
interference pattern [3l| and is a source of entanglement 
of the continuum wave functions [32]. 

The Hermitian mixing of states tends to make decay 
widths more similar (width equilibration), whereas the 
anti-Hermitian mixing repels widths leading to the for- 
mation of different time scales in the system (width bifur- 
cation). In particular, the anti-Hermitian mixing tends 
to concentrate the continuum coupHng in a single state 
(aligned state) of the OQS. As discussed below, proper- 
ties of this state are crucial for the appearance of cluster 
correlations close to the cluster-decay threshold. 

In the following, we shall employ SMEC to study the 
collectivization of SM states around the threshold. The 
detailed description of SMEC can be found elsewhere [si 
[33| (see also Ref. [l6|). The Hilbert space is divided 
into orthogonal subspaces where index /i denotes the 
number of particles in the scattering continuum {ja = 
0, 1, . . . ). An OQS description of internal dynamics in Qo 
includes couplings to the environment of decay channels, 
and is modelled by the energy-dependent Hamiltonian: 

n{E) = Ho^ Hi{E) =Ho^ VSh{E). (1) 

In this expression, Hq is the CQS Hamiltonian (the SM 
Hamiltonian), Vq is the continuum-coupling strength, E 
is the scattering energy, and h{E) is the coupling term be- 
tween localized states (Qo) and the environment of decay 
channels (Qi, Q2, • • • )• The 'external' mixing of two SM 
eigenstates i and j due to Hi{E) consists of the Hermi- 
tian principal value integral describing virtual excitations 
to the continuum and the anti-Hermitian residuum that 
represents the irreversible decay out of the internal space 
Qo- The SMEC solutions in Qo are found by solving the 
eigenproblem for the non-Hermitian Hamiltonian ([1]), 

^QoQol^.) = ^,(^,^0)1^,) 

(^l-^QoQo = E;{E,Vo){^jI (2) 

in the biorthogonal basis {^j\^k) = ^jk- The left 
and right (<I>j| eigenstates are related by complex conju- 
gation. 




FIG. 2: The continuum-coupling energy correction JSj) to the 
SM eigenstates as a function of energy above the proton decay 
threshold Ethr- The symbols ii, 12 denote two eigenstates 
within the energy interval characteristic of a strong coupling 
to the common decay threshold. The Coulomb interaction 
suppresses the continuum coupling and shifts the maximum 
of the correlation energy, Etp, above the decay threshold. 

An explicit energy dependence of the effective Hamil- 
tonian ([1]) is the origin of its strong non-linearity. More- 
over, the continuum-coupling term Hi{E) generates ef- 
fective many-body interactions in the internal space, even 
if it has two-body character in the full space. For in- 
stance, the continuum-coupling correction to the bind- 
ing energy induces the change in the effective two-body 
monopole terms that cannot be distinguished from the 
effect of a three-body force in SM calculations [34]. 

The eigenstates of H are mixtures of SM eigen- 
states 

l^.)^l*i) = E-^^-'IV'i)' (3) 

where hji is an orthogonal transformation matrix. The 
squared matrix element is the (complex) weight of the 
SM eigenstate ifji in the correlated OQS eigenstate |<I>j): 

Y.TZ{h%) = 1, ^ = 0. (4) 



A. Continuum-coupling correlation energy 

In general, the continuum coupling lowers the binding 
energies of eigenstates that are close to the decay thresh- 
old. The continuum-coupling correlation energy of the 
SM eigenstate {ipi)^ 

eH-AE) = mn - Horn ^ Vo^{^Pi\h{E)m, (5) 

depends on the structure of the SM eigenstate and the 
nature of the decay channel. The correlation energy 
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is peaked at the threshold only for the coupling to the 
i = neutron decay channel [3 5]. For higher partial waves 
and/or for charged particles, the centrifugal Coulomb 
barriers shift the maximum of the correlation energy 
above the threshold. 

Figure [2] shows a typical shape of the continuum- 
coupling correlation energy (j5j) for SM states that are 
coupled to the proton decay channel. The correlation 
energy is proportional to Vq and to the average density 
of SM states in the energy interval £:/2 on both sides of 
the decay threshold. An interplay between the Coulomb 
interaction and the continuum coupling determines the 
peak of the continuum coupling correlation energy, i.e., 
the full width at half maximum a and the centroid E^p of 
the correlation energy around the maximum. The maxi- 
mal value of the continuum coupling energy at the max- 
imum can reach several MeV, depending on the configu- 
ration of SM states involved and the nature of the decay 
channel. Typically, a of the correlation energy peak is 
less than ~0.5 MeV in case of the coupling to the charge 
particle decay channel. This is much less than the width 
of the correlation energy for the coupling to the neutron 
decay channel, which varies in the range 1 MeV< a <3 
MeV [33, [111 . The narrowing of a in the charged particle 
case is a consequence of the Coulomb barrier that cuts off 
the low-energy tail of the continuum-coupling correlation 
energy. 

Even though the continuum-coupling correlation en- 
ergy is a small fraction of the total binding energy, it is 
of the same order of magnitude as the pairing correla- 
tion energy, which profoundly impacts the independent 
particle motion in the atomic nucleus. It is therefore 
pertinent to ask whether correlations induced by the con- 
tinuum coupling may lead to instabilities of certain SM 
eigenstates near the channel threshold, and whether such 
couplings are capable of generating collective rearrange- 
ments involving many SM eigenstates. This question will 
be addressed in Sees. HVl and IVl below. 

An illustration of the continuum-coupling-induced col- 
lective features of SM states is the distribution of weights 
of SM states for a given OQS eigenstate. Figure [3] shows 
real parts of weights b'j^ for the two 0^ SMEC eigenstates 
of ^^Ne. In this calculation, we consider four = 0^ 
SM eigenstates coupled to the i = proton decay channel 
leading to the 1/2+ state of ^^F. The SMEC calculations 
have been carried out in the (Opi/2, 0(i5/2, lsi/2) model 
space, in which there are only four 0+ states in ^^Ne. 
For Hq we take the ZBM Hamiltonian [36] . The residual 
coupling between Qo and Qi is generated by the contact 
force V12 = VoS{ri — r2). Since the reference energy scale 
in SM is undefined, we put the zero energy in SMEC at 
the lowest particle emission threshold. 

The size of the continuum-coupling correlation energy 
depends on the continuum coupling strength Vq , the den- 
sity of CQS states, and the system excitation energy E. 
The range of relevant Vq values can be determined, e.g., 
by fitting decay widths of the lowest states in ^^F. For a 
ZBM Hamiltonian, the experimental decay widths of the 
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FIG. 3: The real part of the weights b% gj for 0+ (top) and 0+ 
(bottom) SMEC eigenstates in ^^Ne. In both cases, the OQS 
eigenstates are at the maximum of the continuum-coupling 
correlation energy ^tp- The continuum-coupling strength is 
Vo = -IGeV fm^ (left) and Vo = -3GeV fm^ (right). 



ground state 1/2^ and the first excited state 5/2^ in -^^F 
are reproduced by taking Vq = —3500 ± 450 MeV-fm^ 
and Vb = -1100 ± 50 MeV-fm^, respectively. The error 
bars in Vq reflect experimental uncertainties. 

The weights lZ{h\^) are shown in Figs.lSJa) and[3](b) for 
Vo = — IGeV fm^ and Vq = — 3GeV fm^, respectively. In 
this calculation, the lowest eigenstate $1 has an energy 
E = £^TP- The largest weight in <I>i corresponds to the 
lowest energy SM eigenstate The weights of higher 
SM eigenstates rapidly drop with i. The distribution 
Tl{h\^) becomes broader if the strength of the continuum 
coupling, or equivalently the density of SM eigenstates, 
is enhanced. The width of this distribution is a measure 
of the collective response to the continuum coupling. 

The distribution of weights 7^(&2i) ^^r the second eigen- 
state ^2 is plotted in Figs.[3Kc) andlS^d). In this case, the 
eigenstate ^2 at £^ = £^tp is close to the threshold and 
<I>i is well bound. Here, the largest weight corresponds 
to ?/^2- With increasing continuum-coupling strength, the 
weights of 7^(^22) ^(^21) become comparable. 



B. Double poles of the scattering matrix 

For E <{) (bound system), eigenvalues Ei{E) of the ef- 
fective Hamiltonian 1-L{E) ([T]) are real. In the continuum 
region, Ei{E) correspond to the poles of the scattering 
matrix and becomes complex-symmetric. The com- 
petition between Hermitian and anti-Hermitian parts of 
1-L{E) may lead to the coalescence of two eigenvalues, i.e., 
to the formation of EPs (double poles of the scattering 
matrix) [28]. The relation of EPs to avoided level cross- 
ings and spectral properties of Hermitian svstemsjs?!, Issf 
has been discussed in considerable detail in schematic 
models [39]. 

EPs are indicated by common roots of the two equa- 
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FIG. 4: Different phases in the scattering continuum of ^Be. 
Of particular interest are two doublets of nearly degenerate 
2^ and 3^ resonances in the vicinity of ^Li+p and ^Be+n 
thresholds, respectively. 



tions: 

— det [n {E; Vo) -SI]=0, u = 0, 1. (6) 

Single-root solutions of Eq. (|6|) correspond to EPs asso- 
ciated with decaying (capturing) states. For the Hamil- 
tonian ([1]), the maximum number of such roots is 

Mmax = 2n{n-1), (7) 

where n is the number of states of given angular momen- 
tum J and parity tt. The factor 2 in Eq. ([7]) comes from 
the quadratic dependence of H on ([1]). 

An excellent example is the 2+ doublet in ^Be 
shown in Fig . 21 which was extensively discussed in 
Refs. |o|, Ef. These two resonances, which have ap- 
proximately the structure of ^Li+p and ^Be+n threshold 
configurations [4^, are seen below one-nucleon emission 
thresholds. Above these thresholds, one finds a nearly 
degenerate doublet of 3+ resonances. 

The nucleus ^Be is an example of how the presence of 
the reaction threshold can infiuence the nature of regular 
solutions. Above the a -\- a threshold, the ground state 
0+ and rotational excitations 2+, 4+ all belong to the 
the same analytic phase which carries the imprint of the 
alpha-decay channel. The new phases start slightly be- 
low the ^Li+p threshold and continue above the ^Be+n 
threshold. The isospin symmetry of the nuclear force and 
the approximate symmetry of these two configurations 
with respect to the Coulomb force determine the strong 
mixing of resonances in these two phases and give rise 
to the formation of energy doublets that resemble ^Li+p 



and ^Be+n mirror thresholds. These 2+ and 3+ doublets 
are thus convincing examples of avoided crossing of res- 
onances, the mechanism of which has been explained in 
a simple two-resonance picture [io]. 

The mixing of wave functions in CQSs is closely re- 
lated to the features of EPs of the complex-extended 
CQS Hamiltonian. Similarly, to understand the mixing 
of wave functions in OQSs one should analyze the spec- 
trum of EPs of the complex-extended OQS Hamiltonian. 
Since the OQS Hamiltonian ([1]) is energy dependent, the 
essential information about configuration mixing is con- 
tained in exceptional threads, ETs, the trajectories of 
coalescing eigenvalues Ei^{E) = Ej^jE) of the effective 
Hamiltonian for a complex value of Vq [43| . 



IV. CONFIGURATION MIXING NEAR THE 
CHARGE PARTICLE EMISSION THRESHOLD 

A complete picture of the near-threshold configuration 
mixing in CSM many-body wave functions can be ob- 
tained by investigating properties of ETs of the complex- 
extended CSM Hamiltonian (2H . The anti-Hermitian 
mixing leads to an accumulation of collective effects in a 
single state close to the particle-emission threshold. Such 
an aligned state couples strongly to the decay channel 
and, therefore, carries many of its characteristics. A large 
part of the total continuum-coupling correlation energy 
is concentrated in the aligned state [21| as can be seen in 
the energy dependence of ETs in Fig. [5l 

Figure [5] displays ETs for the 0+ eigenvalues of the 
complex-extended SMEC Hamiltonian in ^^Ne as a func- 
tion of |Vb|. There are four 0+ SM eigenstates in this 
model space and, hence, six different ETs for decaying 
resonances. All ETs which are relevant for the collec- 
tive mixing of SM states involve the aligned eigenstate 
Ot- These ETs: 0]^ - Oj, 0]^ - Oj, and - Oj are 
shown in the lower panel. They exhibit a minimum of 
I Vol (the turning point, TP) at the same proton energy 
Etp — 0.55 MeV. At this energy, the collective mixing of 
SM states is maximal. The remaining ETs, which do not 
involve the aligned state, are found at exceedingly large 
values of I Vol (see the upper panel in Fig. [5]) and have no 
infiuence on the configuration mixing. 

We assert that the collective mixing of all SM eigen- 
states with the same quantum numbers via the 
aligned state is the essence of the clustering phenomenon. 
This mixing, which results from the opening of the many- 
body system, can explain the appearance of cluster states 
close to the corresponding cluster decay threshold. The 
aligned state is also an archetype of the cluster state. For 
charged-cluster configurations, the appearance of cluster- 
ing depends on whether an SM state is found inside of 
the energy window given by the competition between the 
continuum coupling and the Coulomb interaction (see 
discussion around Fig. [2]). Consequently, the window 
of opportunity for charged-particle clustering is situated 
around the turning point above the charged-cluster decay 
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FIG. 5: ETs for the 0^ eigenvalues of the complex-extended 
effective Hamiltonian ([T]) in ^^Ne plotted as a function of |Vb|. 
Only ETs corresponding to decaying resonances are shown. 
Each point along an ET is an EP labelled by the proton energy 
with respect to the ^^F(l/2+)+p(-£ = 0) threshold. There 
are six threads in this case. Three of them (bottom panel), 
involving the aligned 0^ state, are found in the physical range 
of I Vol values and correspond to a coalescence of 0^ — J (solid 
line), 0^ — 0^ (long-dashed line), and 0^ — (short-dashed 
line) eigenvalues. The remaining three threads (top panel), 
~ — O4 , and 0^ — O4 , correspond to extremely 

large values of |Vo|. For more details of the model space and 
interaction, see Sec. IIIIl 



threshold. 



V. MAXIMUM CONTINUUM COUPLING 

The mixing of SM eigenstates through the continuum 
is strongest in the narrow window of proton energies 
around the turning point for all ETs which involve the 
aligned eigenstate. In the following, we shall illustrate 
some pertinent properties of the turning point. To this 
end, we consider 0+ SM states in ^^S, which are cou- 
pled to the same one-proton decay channel ^^P(l/2^) 
+ p{i = 0). SMEC calculations in ^^S have been per- 
formed in the (151/2,0^5/2,0(^3/2) SM model space us- 
ing the USDB interaction [4^. Since we consider two 0^ 
states, there are two ETs associated with decaying reso- 
nances. 

Figure [6] shows the energy and width of two 0+ eigen- 



> 







-0.5 




0.5 0.2 0.4 0.6 

proton energy [MeV] 



FIG. 6: Energies Ei (left) and widths F^ (right) of 0+, {i = 
1,2) eigenvalues of the effective Hamiltonian ([T]) in ^^S as a 
function of the proton energy in the vicinity of the ^'^P(l/2^) 
-h p(-£ = 0) decay threshold. The aligned state is shown by a 
solid line. Calculations have been performed using the SMEC 
Hamiltonian with Vb = — 1 GeV fm^. The threshold energy 
Ethr is put at E = 0. 



values of V, in the vicinity of the 2^P(l/2+)+p(^ = 
0) decay threshold. The continuum coupling has the 
strongest effect on the lowest-energy eigenvalue Ei . This 
aligned eigenstate becomes fairly broad above the pro- 
ton emission threshold. There appears a narrow min- 
imum in Ei{E) around the turning point energy ^tp- 
In general, Etp depends on the height of the effective 
(Coulomb+centrifugal) barrier, the structure of SM wave 
functions involved, and the decay channel itself. 




-1 1 

proton energy [MeV] 

FIG. 7: Continuum-coupling correlation energy for the 
aligned 0^ state in ^'^S as a function of the proton scat- 
tering energy for various Coulomb potentials. The solid 
Une corresponds to the Coulomb potential of ^^P+p system 
{Zt = Zi • Z2 = 15; Z = Zi + Z2 = 16) and a uniformly 
charged sphere of radius i?o = 4 fm. The dashed lines corre- 
spond to artificially reduced Coulomb potentials with Zt = 8 
and radii Ro = 3.61 fm (short-dashed line) and 4 fm (long- 
dashed line). The latter radius is calculated from the rela- 
tion Ro = 1.35 (A — 4)^/^ + 1.3 fm and represents the ra- 
dius of the ^Be + a potential (Zt = 8; Z = 6). Calcula- 
tions have been performed using the SMEC Hamiltonian with 
Vo = -1 GeV fm^ 

Figure [3 shows the continuum-coupling correlation en- 
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ergy ([5]) for the aligned 0+ SMEC eigenstate in ^^S, which 
is calculated for various Coulomb potentials. The solid 
line is calculated for the Coulomb potential of ^^P+p 
system {Zt = Zi • Z2 = 15; Z = Zi -\- Z2 = 16, where 
Zi = 1, Z2 = 15 are proton numbers of the colliding 
nuclei.). Other curves have been calculated using the 
same continuum-coupling matrix elements of ^^S and the 
arbitrarily reduced strength of the Coulomb potential 
{Zt = 8) to assess the impact of the Coulomb potential 
on the continuum-coupling correlation energy. It is seen 
that the correlation energy decreases rapidly with the in- 
creasing Coulomb barrier. For the same value of Zf of 
fragment charges, the magnitude of the correlation en- 
ergy decreases with the increasing radius of the potential 
Rq. Even though the width of the energy window cor- 
responding to large continuum-coupling correlation en- 
ergy increases with increasing Z^, it remains significantly 
smaller than in the neutron case [35]. 

Figure [8] shows the energy dependence of the off- 
diagonal matrix element 612 (|3]) for the 0+ aligned state 
in ^^S. The minimum of bi2{E) corresponds to the pro- 
ton energy at which the largest mixing takes place. This 
energy is nearly equal to Etp- Energy variations of 612 
with Z and Rq are analogous to those found for £^corr (cf. 
Fig. [7]). Notice that the magnitude of mixing, as given 



-0.005 




-0.004 I ^ ' ^ J — ^ ' 1 

-10 12 

proton energy [MeV] 

FIG. 8: Real and imaginary parts of the amplitudes bji of 
the orthogonal transformation Q for the aligned 0+ state 
in ^^S plotted as a function of the proton scattering energy 
for different Coulomb potentials as in Fig. [71 The point of 
maximum mixing at Etp is marked by arrows. 

by 1 612 I, is strongly reduced with increasing proton 



number Z. Consequently, the imprint of the branch- 
ing point on regular solutions of the scattering problem 
involving charged particles is reduced in heavier nuclei. 
This generic feature, with far-reaching consequences for 
the manifestation of nuclear clustering, does not depend 
on the particular emission channel and the angular mo- 
mentum i involved. 

The behavior of Ecott{E) and bi2{E) around the turn- 
ing point energy Etp follows from properties of ETs 
of the complex-extended CSM Hamiltonian. Figure [9] 
shows the ETs in the complex- plane calculated for 
the aligned 0+ state in ^^S assuming different Coulomb 
barriers. For each ET, one can see two regions where the 
trajectory turns. The first one, at large values of |7^(Fo)| 
and is seen at the threshold energy E = Efhr- 

The second one, at smaller values of Vq^ is found at the 
energy E = E^p. It is obvious from the X{Vo) — 1Z{Vo) 




-400 -200 

9t Vo [MeV-fm^] 



FIG. 9: ETs in the complex- Vb plane for the aligned 0^ state 
in ^^S. SMEC calculations have been performed for different 
Coulomb barriers as indicated. The solid line shows results 
for the standard Coulomb barrier; in the variants marked by 
dashed lines the Coulomb barrier has been reduced. 

plot that the continuum-coupling induced mixing of SM 
states diminishes with increasing Zt because the turning 
point at ^ = ^TP moves farther away from the physical 
limit (XVb=0) towards larger values of |Vb|. This generic 
tendency does not depend on the nature of the decay 
threshold. In other words, the increase of the Coulomb 
barrier in heavier nuclei not only suppresses correlation 
energy close to proton-decay threshold, but also prevents 
other kinds of clustering such as the alpha-clustering near 
the alpha-emission threshold. Conventionally, the lat- 
ter effect is associated with the disruptive influence of 
the nuclear spin-orbit interaction on alpha-clustering (see 
Ref.j45| and references cited therein), i.e., with the spe- 
cific feature of the nuclear potential. 

Whereas, for a fixed value of Zt, the energy Etp of the 
turning point is independent of the continuum-coupling 
constant Vb, the energy E^ax corresponding to the max- 
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imum continuum correlation energy depends on Vq that 
is adjusted to the experimental data on low- lying reso- 
nances. The mismatch between these two energies in- 
creases when the 'physical' continuum-coupling strength 
becomes significantly different from |1^|(£^tp), i.e., at 
large values of Zt. 

Since Vq determines the overall size of the continuum- 
coupling term, the optimal energy window for the ap- 
pearance of cluster correlations should be largely inde- 
pendent of the detailed structure of Hi{E) ([1]). In other 
words, given a value of Z^, one expects the turning point 
energy ^tp to be weakly dependent on the nature of the 
charged particle decay channel and the parameters of the 
potential. Only the states in the vicinity of the turning 
point at Ftp — Emax may exhibit cluster correlations. 
Whether they actually do or not depends mainly on a 
few non-generic parameters such as Vq or the radius of 
the interaction potential. 




5 10 15 

z, 

FIG. 10: Turning point energy ^tp (bottom) and the ratio 
of Ftp to the height of the Coulomb barrier Eb (top) versus 
Zt for the aligned 0^ states in ^^S (lines) and ^^Ne (crosses). 
SMEC calculations are performed for Rq = 3.61 fm (solid line) 
and 4fm (dashed line). 

To see the impact of the Coulomb barrier and the po- 
tential radius on £^tp, the lower panel of Fig. [TOl shows 
the dependence of E^p on Zt in ^^S calculated in SMEC 
for two values of Rq and varying arbitrarily the strength 
of the Coulomb potential. The point Zt = 8 and Rq = 4 
fm corresponds to the ^Be + a case, i.e., is relevant to 
the Hoyle resonance in ^^C. Figure [TOl also indicates ^tp 
for the aligned state 0^ in ^^Ne, which is calculated us- 
ing Ro = 4fm and the reduced charge Zt = 8. This 



point agrees well with the curve E^piZt) obtained using 
^^S wave functions, i.e., J^tp seems to be nearly indepen- 
dent of the potential radius and the detailed structure of 
SM wave functions involved in the continuum coupling. 
The upper panel of Fig. [10] displays the ratio of Etp 
and the height of the Coulomb barrier, Eb^ as a func- 
tion of Zt. Also in this representation, the dependence 
on the potential radius and the detailed structure of SM 
wave functions are rather weak. Both observations sup- 
port our assertion that one may use the curve ETp{Zt) of 
Fig. [To] to read out the position of the centroid for opti- 
mal continuum-induced cluster correlations for any kind 
of charged clusters. 



VI. CONCLUSIONS 

Nuclear clustering is a key problem in low-energy nu- 
clear structure studies. The challenging problem is how 
to separate generic and specific features of this intricate 
many-body phenomenon. The energetic order of parti- 
cle emission thresholds, and their nature, is given by the 
N- and Z-dependence of the nuclear binding energy, i.e., 
it depends on specific properties of the nuclear Hamilto- 
nian. The appearance of specific decay channels involving 
both kinds of nucleons, and the absence of stable clus- 
ters/nuclei entirely composed of like nucleons, is a direct 
consequence of the isospin structure of the nuclear force. 
The continuum coupling effects play a rather minor role 
in the interplay between various contributions to the to- 
tal binding energy, at least close to the valley of stability. 

On the other hand, the phenomenological rule that 
cluster correlations are seen in the vicinity of the re- 
spective cluster emission threshold is unlikely a conse- 
quence of specific properties of nuclear forces and calls 
for a generic explanation. A plausible explanation of this 
rule, put forward in Ref. [2l|, has been in terms of the 
collective coupling of SM states via the decay channel(s). 
This anti-Hermitian coupling leads to the formation of 
the aligned state, the eigenstate of the OQS which cap- 
tures most of the continuum coupling, and, above the 
decay threshold, exhausts most of the decay width. 

The mechanism responsible for the creation of an 
aligned near-threshold state is mathematically similar to 
the formation mechanism of collective super-radiant or 
trapped states [T2|-[l4j. However, the physical domain of 
aligned states is not restricted to the region of large den- 
sity of resonances in the continuum. It can even corre- 
spond to a bound state at energy below the lowest decay 
threshold. 

The collectivity of an aligned state is a fingerprint of 
instability in an ensemble of all SM states having the 
same quantum numbers and coupled to the same decay 
channel [23. Quantitatively, manifestations of this insta- 
bility depend on the strength of the continuum coupling, 
the density of SM states, and the nature of the decay 
channel. For large angular momenta i and/or charged 
particle decays, the energy window for the formation of 
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a collective aligned state is pushed by centrifugal and 
Coulomb barriers above the decay threshold. 

Are the examples presented in this study, e.g., = 
0+ wave functions of ^^^S affected by the 23p(i/2+) + 
p(^ = 0) proton decay channel, relevant to near-threshold 
clustering effects involving composite clusters? We have 
demonstrated that the point of maximum continuum cou- 
pling for the charge-particle decay is related to the turn- 
ing point of all ETs involving the aligned state. The 
energy of this turning point is predicted to be fairly in- 
dependent on both the continuum-coupling strength and 
the detailed structure of coupling matrix elements but it 
varies with the potential radius and the product of frag- 
ment charges Zt. This means that one can in principle 
estimate the position of the energy interval of maximum 
continuum coupling for any charge-particle decay channel 
using the above results for the one-proton decay. On the 
contrary, the quantitative description of decay properties 
of the cluster states is more involved and requires a mi- 
croscopic calculation of both the coupling to the cluster 
decay channel(s) and the cluster formation amplitude. 

The continuum-coupling correlation energy and the 
collectivity of mixing amplitudes are reduced with in- 
creasing Z and the charge of emitted particle. This re- 



sults in a narrow energy window for the appearance of 
collective alpha-cluster correlations in medium and heavy 
nuclei. 

In the studied case of one-proton decay, the aligned 
state becomes a broad resonance already at the border 
of the 5(i-shell. This problem is less acute for alpha decay 
because of the small preformation factor in question [25l. 
46]. Consequently, very narrow aligned alpha- decaying 
states can be found around the turning point in quite 
a few p- and sd- shell nuclei. The Hoyle resonance in 
^^C represents a splendid example of such a continuum- 
correlated collective state. 
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